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STATIC- TEMPERATURE DISTRIBUTION IN A 


FLAT-PLATE COMPRESSIBLE TURBULENT BOUNDARY LAYER 

WITH HEAT TRANSFER 

By S. Z. Pinckney 
Langley Research Center 

SUMMARY 

An expression based on the differential equations for the local energy transfer and 
shear is derived for the relationship between the static temperature and the velocity in a 
zero-pressure-gradient turbulent boundary layer. The boundary conditions imposed by 
these governing equations showed the temperature-velocity relation to be at least a fourth 
degree polynomial which differs considerably from the second degree polynomial obtained 
by previous investigators. These governing equations are integrated by using as bound- 
ary conditions the value of a particular total-energy-deficiency parameter of the boundary 
layer at the station under consideration, the slope of the laminar- sublayer temperature- 
velocity profile at the wall (as determined from a modified Reynolds analogy), and the 
local free-stream and wall conditions. Results obtained by this method are compared 
with experimental temperature-velocity profile data and are found to correlate well with 
the experimental profiles. 


INTRODUCTION 

The ability to predict the variation of turbulent-boundary-layer parameters, such 
as momentum, displacement, and energy thicknesses, is very important in supersonic 
flow. In order to describe the distribution of these parameters, a knowledge of the 
turbulent-boundary-layer temperature profile is needed, in addition to the velocity profile 
and friction coefficients. 

The difficulties encountered in the theoretical evaluation, even for the zero- 
pressure- gradient turbulent boundary layer, have resulted in the extensive use of empiri- 
cal correlations. For example, correlations for the compressible friction coefficient 
(ref. 1) and zero-pressure-gradient incompressible velocity profiles (ref. 2) have 
resulted in various empirical methods for the determination of these parameters. 

Up to the present time, the temperature profile for the compressible turbulent 
boundary layer with heat transfer has not been correlated accurately even for zero 



pressure gradient. The determination of temperature profiles has been limited, gener- 
ally to attempts at altering an expression obtained by Luigi Crocco for laminar flow (this 
expression is presented in ref. 3) . In reference 4, however, an expression for the tem- 
perature profile that is based on the differential equations for local energy transfer and 
shear is presented. 

In the determination of the expression for the turbulent-boundary- layer tempera- 
ture profile, it was assumed in reference 4 that the ratio of local energy transfer to wall 
heat transfer divided by the ratio of local shear to wall shear is a constant. Also the 
ratio of eddy diffusivity of heat to eddy diffusivity of momentum was assumed to be a 
constant across the boundary layer. It is pointed out in reference 4 that for flow over a 
flat plate this ratio of eddy diffusivities is approximately unity; therefore, most of the 
calculated results presented in reference 4 were derived with a value of 1 assumed for 
the ratio. In the present paper, a comparison of the results obtained by using the theo- 
retical expression for the temperature profile of reference 4 with experimental data 
(ref. 5) reveals substantial differences between the theory and some flat-plate tempera- 
ture profiles. This observed discrepancy between theory and experiment is believed to 
be due in part to the aforementioned assumptions, as well as to other simplifying assump- 
tions. The method developed in the present paper avoids the difficulties of evaluating 
local values of some of the less tractable boundary- layer parameters by making use of 
an integral parameter which exerts a strong influence on the temperature profile. 

The present analysis is similar in approach to that presented in reference 4, but a 
polynomial expression in terms of the ratio of local velocity to free-stream velocity is 
assumed for the ratio of local energy transfer to local shear, and a function which con- 
tains eddy diffusivities of heat and momentum as parameters is determined from certain 
boundary conditions. In particular, a boundary condition is imposed that requires the 
value of a particular total-energy-deficiency parameter of the boundary layer, deter- 
mined from the theoretical static-temperature — velocity profile relation, be equal to the 
known value of the total-energy-deficiency parameter of the boundary layer at the station 
under consideration. The use of a total-energy-deficiency parameter is a new and 
reasonable approach to the problem based on the observation of many researchers that 
most turbulent-boundary-layer parameters are functions of the previous history of the 
boundary layer. At the present time, for example, in a design problem the lengthwise 
distribution of the total energy deficiency and the momentum thickness needed to deter- 
mine the value of the total-energy-deficiency parameter could be determined by integra- 
tion of the heat transfer to the wall and of momentum losses due to shear along the wall, 
respectively. In a converse sense, if the total energy deficiency and momentum thickness 
cannot be determined at a given station, the comparisons of experimental data and results 
from theories presented herein strongly imply that within the current state-of-the-art 
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the static-temperature — velocity profile cannot be determined with accuracy except for 
certain cases. 
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SYMBOLS 


7T- + 


expression for 


pe h 


l± + pe 

constants assumed in the relationship between q/r and u/u 0 

2t w 

friction coefficient, tj- 

P oo U - 

specific heat of a fluid at constant pressure 
static enthalpy (for perfect gas, h = CpT) 
thermal conductivity 
free-stream Mach number 

exponent for y/ 6 in velocity profile relation (^- = ^ ^ 


Prandtl number, 


c p P 


energy transfer in y-direction, per unit time per unit area (as defined by 
eq. (2)) 

total energy deficiency of boundary layer relative to free stream, 

,6 


Q= I 0 PU (\oo - h t) dy 


perfect gas constant 
static temperature 
total temperature 
hypothetical wall temperature 
velocity component parallel to plate 
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X 


longitudinal distance along plate 


y perpendicular distance from plate 

y ratio of specific heats, taken as 1.4 for air 

6 boundary-layer thickness 

e eddy diffusivity of momentum 

e ^ eddy diffusivity of heat 

C ® pu / 

6 momentum thickness, \ (1 

^0 Pco U ° oV 

ju viscosity 

p density 

r shear stress, force per unit area 

Subscripts: 

aw adiabatic wall conditions 

t total conditions 

w wall conditions 

°° local free- stream conditions 

ANALYSIS 

Classically, the turbulent boundary layer is divided into three regions: 

(1) The inner region or laminar sublayer; 

(2) The transition region between the laminar sublayer and the turbulent flow 
region; and 

(3) The turbulent flow region. 
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In deriving an expression for the turbulent-boundary-layer enthalpy-velocity dis- 
tribution it was not necessary to distinguish between these three regions of the boundary 
layer. However, various features of two of these regions, namely the laminar sublayer 
and the turbulent flow region, are in evidence through the boundary conditions and the 
form of the equations for shear stress and energy transfer that are assumed. 

From reference 4 the local shear and energy transfer in the boundary layer can be 
expressed as follows: 

r = (m + pe)% (1) 


* = -^ + '* h )ar < “ tP£>u ! (2 > 

Dividing equation (2) by equation (1) and nondimensionalizing h and u with respect to 
the free stream give 

. £- + pe h d (r L ') 

\JW _u_ / 3 n 

r Uoo n + pe d ru_\ 00 Uoo v ’ 

W 


In order to integrate equation (3) a knowledge of the relationships between q/ r and u 

^r- + pe h + P e h 

and between P — and h is needed. Herein the ratio -2 is assumed to be 

p + pe p + pe 

constant and is designated a. 

The following third degree polynomial in terms of u/uoo is assumed for the rela- 
tionship between the ratio q / t and the velocity profile: 




From a consideration of the boundary conditions imposed on the governing equa- 
tions (which are presented in appendix A of ref. 4) when a constant wall temperature is 
assumed, , . . . 

fig! -&) -0 <5a> 


8 y/w Vfy 


Applying the boundary conditions of equations (5a) and (5b) to equation (4) gives d = 0 
and e = 0. Therefore, the polynomial assumed for the relationship between q/ r and 
the velocity profile reduces to 
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( 6 ) 


q 

T 



When a is substituted for 


h + pe h 


and the expression for q/ r from equa- 

li + pe 

tion (6) is used in equation (3), integration from the wall to an arbitrary point in the 
boundary layer gives, after rearranging, 
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(7) 


The coefficients a, b, and c in equation (7) are constant at a given station but vary 
with x; the evaluation of the constants is accomplished by using the following boundary 
conditions: 

(1) From Colburn's version of Reynolds analogy, 

_ h aw ~ Nv Pr l/3 

hoo 

(2) At the outer edge of the boundary layer, — =1.0 and — =1.0; and 

hoo U oo 

(3) The integral of the total energy deficiency across the boundary layer evaluated 
by combining the enthalpy relation derived from equation (7) with a suitable assumed 
velocity profile must equal the total enthalpy deficiency of the entire boundary layer up 

to the x-station under consideration. A derivation for the total-energy-deficiency param- 
eter is presented in appendix A. 

The differentiation of equation (7) and the evaluation of the result at the wall give 



The application 
tion (8) into the 
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of boundary condition (2) to equation (7) and the substitution of equa- 
result give 
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The substitution of equations (8) and (9) into equation (7) leads to 
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and the rearrangement of equation (10) using boundary condition (1) gives 
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With the assumption of a perfect gas (for which h = CpT^, equation (11) can be 
rearranged to give 


(12) 


In equations (11) and (12) the only remaining unknown is the ratio of constants c/a. 
The final form of equation (11) or (12) and the value for the ratio c/a are determined 
through application of boundary condition (3). This boundary condition requires the 
matching of the value of the total- energy-deficiency parameter determined by using equa- 
tion (A3) in conjunction with the theoretical enthalpy-velocity or temperature-velocity 
profile relation (eq. (11) or (12), respectively) with the known value of the total-energy- 
deficiency parameter of the boundary layer at the station under consideration. This 
matching process and thus the determination of the correct value for the ratio of con- 
stants c/a consist of the substitution of the theoretical enthalpy-velocity or 
temperature-velocity profile relation along with a suitable velocity profile into the right- 
hand side of equation (A3) and iterating, with several values tried for the unknown c/a. 
This iteration is complete and the correct value for the ratio of constants c/a is 
determined when the right-hand side of equation (A3) produces a value for the total- 
energy-deficiency parameter equal to the known value. In practical application of the 
present theoretical method, the total energy deficiency and the momentum thickness 
needed to determine the known value of the total-energy-deficiency parameter would be 
determined from the heat loss to the wall and friction losses along the wall integrated 
from the boundary- layer origin to the x- station under consideration. However, in the 
present comparisons of experimental data and theoretical results, sufficient informa- 
tion for performing this type of integration accurately generally was not available; there- 
fore data were chosen for which the total temperatures across the boundary layer were 
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available (see table I) and the known value of the total- energy-deficiency parameter could 
be determined from the experimental profiles. 

RESULTS AND DISCUSSION 


Comparisons of Theory and Experiment 


The present theoretical results have been compared with 82 experimental 
temperature- velocity profiles obtained from references 5 to 14. The results of these 
comparisons are presented in figures 1 to 9. For these comparisons a power-law veloc- 


ity profile of the form — = (t) was assumed. Theoretical temperature- velocity 
profiles were obtained for comparison with each experimental temperature-velocity pro- 
file by using a high-speed digital- computer program, with values of 5, 7, and 9 assumed 
for the power-law velocity-profile exponent N. As exemplified by figures 1(a) and 1(b), 
the theoretical temperature-velocity profile is a weak function of the power-law velocity 
profile assumed. As a result it was considered sufficient to present in figures 2 to 9 only 
the theoretical curves corresponding to a power-law exponent of 5, with the exception of 
figure 8(d). Because the experimental velocity profile for the experimental data of fig- 
ure 8(d) corresponded to an N value of approximately 20, the theoretical temperature- 

u /y 0-/20 

velocity profile of figure 8(d) was also determined with — — = ( - as the assumed 

velocity profile relation. It is evident that the 1/5-power profile is not as adequate an 
approximation for these experimental data as it was for the other experimental data 
presented. 


1/N 


For comparison purposes figures 1(a) and 1(b) include sample theoretical 
temperature-velocity profiles calculated by using the theoretical method of reference 4 
as well as sample profiles calculated by using the modified Crocco distribution assumed 
in reference 15. Figure 10 shows sample experimental velocity profiles which corre- 
spond to some of the experimental temperature -velocity profiles; also shown in figure 10 
for comparison purposes is a 1/7-power velocity profile. 


The experimental temperature-velocity profiles considered in this presentation are 
for a range of M m from 0.85 to 8.18 and a range of T w /Too from 0.997 to 7.38. In 
general, comparisons between results of the present theory and experimental data showed 
maximum discrepancies between point values of temperature of approximately 5 percent 
for the outer 90 percent of the boundary layer (values of y/5 between 0.1 and 1.0) and 
approximately 10 percent for the inner 10 percent of the boundary layer (values of y/ 5 
less than 0.1). In comparisons involving the data of reference 8 (fig. 4) the discrepancies 
are larger, approximately 19 percent at distances from the wall of 0.3 millimeter - that 
is, for y/5 = 0.05. The discrepancies are within 10 percent at values of velocity ratio 
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greater than 0.75 or at distances from the wall greater than 0.5 millimeter - that is, for 
y/6 > 0.08. The reason for the larger discrepancies between the results of the present 
theory and the data from this reference is not known; however, it was indicated in refer- 
ences 16 and 17 that the experimental values of Cj in reference 8 are considered unre- 
liable. Reference 17 also indicated that the profile data of reference 8 may have been 
taken in a transition region. In particular, it should be noted that the experimental 
temperature-velocity data of reference 8 seem to show inconsistencies between profiles. 
For example, the experimental temperature-velocity profiles of figures 4(m) and 4(n) have 
essentially the same boundary conditions: 


Figure 

Mco 

Uoo, 

meters/second 

Reynolds number 
based on 9 

Stagnation 

pressure, 

atm* 

T-w/Toc 

4(m) 

5.26 

791.0 

3880 

7.31 

5.463 

4(n) 

5.29 

786.6 

4300 

7.09 

5.564 


*1 atmosphere = 101 325 N/m^. 


However, a comparison of the experimental total- energy-deficiency profiles corresponding 
to the data in figures 4(m) and 4(n), as presented in figure 11, shows significant differ- 
ences, with the more downstream profile (that corresponding to the data in fig. 4(n)) 
having the lower total energy deficiency. 


In the present comparison of theoretical results and experimental data, the apparent 
absence of any effect resulting from not considering the laminar- sublayer velocity profile 
is believed to be due to the existence of thin laminar sublayers. The thickness of the 
laminar sublayers of the experimental data corresponds approximately to only 5 to 10 per- 
cent of the total boundary-layer thickness and therefore makes little difference to any 
integral across the boundary layer. As a matter of interest, the edge of the laminar sub- 


JL v* r 

layer as estimated by the relation — = 12\— — — (ref. 2) is presented in figure 1. It is 

Uoo \1T oo 2 

believed that the use of an assumed power-law velocity profile in the theoretical calcula- 
tions must be restricted to turbulent boundary layers with relatively thin laminar 
sublayers. 


Index Parameter for Ideal Flat- Plate Flow 

As discussed in the introduction, most previous attempts to express the relation- 
ship for the turbulent-boundary-layer temperature-velocity profile have been restricted 
to modifications of the expression obtained by Crocco for laminar flow. A typical example 
of these temperature- velocity profile relationships (ref. 15) is 
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(13) 



r - i 

2 


Pr^M 2 


1 - 




For certain examples (see fig. 1(b)), experimental data were found to correspond closely 
to results from equation (13) (method of ref. 15) as well as to results obtained by the theo- 
retical method of reference 4. In other examples (see fig. 1(a)), significant differences 
between experimental data and results obtained by using the methods of references 4 
and 15 are evident. It should be noted that the difference in theory and experiment shown 
in figure 1(a) could result in as much as a 50-percent error in a computation of boundary- 
layer displacement thickness. Therefore the theoretical methods of references 4 and 15 
were explored further by examining the derivation of the pertinent governing equations. 
This examination revealed that the theoretical derivations had as a basis the assumption 
of the existence of a constant flat-plate surface temperature and free- stream Mach num- 
ber from the plate leading edge to the station under consideration. Based on these 
assumptions, Colburn's form of Reynolds analogy gives (for a perfect gas) the total energy 
deficiency at the station under consideration as 


Poo Uo oCpToo /T aw - 

Q ■ p r 2/3 V 'T~ / 

or, solved for T^T^, 

T w _ T aw _ Pr 2//3 Q 

Too Too T^PooUooCp^# 


(14) 


(15) 


On the basis of these observations for flat-plate flow, it is believed that good correspon- 
dence between experimental data and the theoretical results of references 4 and 15 can be 
obtained when T w /Too, Q, and 6 are related as indicated by equation (15), which 
corresponds to ideal flat-plate flow. 


From this discussion, it is apparent that an index parameter is needed to provide a 
quantitative value for the deviation of a profile from ideal flat-plate flow. Equation (15) 
may be used for this purpose. For ideal flat-plate flow to exist, the wall-temperature 
ratio T w /Too must equal the difference between the adiabatic-wall-temperature ratio 

Taw/Too and the total-energy-deficiency parameter — (see eq. (15)). The 

' , TooP^ooCp,^ 

Pr 2 / 3 0 

total- energy-deficiency parameter — is directly proportional to the total- 


TooP^ooCp^fi 1 


energy-deficiency parameter 


Q 


used in this report as a boundary condition to 


p u^hj. 6 
^co 00 t,oo 

determine the ratio of constants c/a of equation (11) or (12). When the boundary layer 
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does not correspond to ideal flat-plate flow, the wall temperature T w of equation (15) 
is arbitrarily replaced by the hypothetical wall temperature and thus the following 

equation is generated: 


Tw _ Taw _ Pr^^Q 

Too Too TooP^UooCp^ OQ0 


(16) 


Now T* is defined as the hypothetical wall temperature required under ideal flat-plate 
w 

flow conditions to produce a given value of the total-energy-deficiency parameter. Thus 
the difference between T w and T^ is a measure of the deviation from ideal flat-plate 
flow as well as a function of the total-energy-deficiency parameter. The value of 


L w 


rp ^ 

± W 


along with other pertinent information is presented in table I for each experi- 


w 


mental temperature-velocity profile considered. Since the index parameter 


L w 


w 


rp* 

W 


and the ratio of constants c/a are both functions of the total- energy-deficiency param- 


eter — , a means is suggested for the correlation of the ratio of constants c/ a 

Pco U ooht ? 

needed for use in the present theoretical enthalpy- velocity or temperature-velocity rela- 
tion of equation (11) or (12). Equation (16) was combined with the perfect-gas version of 
the theoretical temperature-velocity relation (eq. (12)) with the use of a high-speed 
digital- computer program, and the set of curves presented in figure 12 were generated. 
The curves in figure 12 are for a power-law exponent of 7 and free-stream Mach numbers 
of 3, 5, and 7. Each curve is for a particular ratio of wall to free-stream temperature. 
Figure 12 can be used to determine a value for the ratio of constants c/a directly, 
instead of by iteration. For example, the integrated total energy deficiency Q and the 
momentum thickness 6 needed for use in equation (16) would be determined from the 
heat loss to the wall and friction losses along the wall integrated from the boundary- layer 
origin to the x- station under consideration. Then use of these values of Q and 6 in 
equation (16) would produce a value of T^. Upon combining the resulting value of T^ 

T - T* 

with the wall temperature T w , the value of — — — needed for entrance into figure 12 

T w 

would be obtained. The resulting value of c/a would be the same as that which would 
be obtained by the iteration process with a l/7-power velocity profile used in the energy- 
deficiency integral of equation (A3). 


CONCLUDING REMARKS 


The present method for computing the static-temperature — velocity distribution for 
a zero-pressure-gradient compressible turbulent boundary layer has been developed by 
using as a boundary condition the known value of a particular total- energy-deficiency 
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parameter of the boundary layer at the station under consideration. In comparing results 
of this analytical method with experimental data, a power-law turbulent-boundary- layer 
velocity profile was assumed. It was observed that the analytical results are a weak 
function of the power-law velocity profile assumed and that the theory produces results 
which agree well with the experimental data. The experimental velocity profiles, for 
which the power-law assumption gives good correspondence, are restricted to turbulent 
boundary layers that have thin laminar sublayers. For boundary layers with relatively 
thick laminar sublayers, more accurate velocity profile representations should be used. 

A comparison of experimental temperature-velocity profiles with modified Crocco 
flat-plate turbulent-boundary- layer temperature-velocity relations revealed that correla- 
tion with experimental profiles is obtained when the experimental data correspond to ideal 
flat-plate flow. The present method gives good agreement for all experimental profiles 
examined, including those for ideal flat-plate flow and those for gross deviations from 
ideal flat-plate flow. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 28, 1968, 

126-15-03-01-23. 
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APPENDIX A 


EQUATION FOR TOTAL ENERGY DEFICIENCY 


The total energy transferred to the wall upstream of a particular boundary- layer 
station can be obtained by the integration of the energy deficiency of the boundary layer 
relative to the free stream. The result is 

p 6 p 6 

Q = J o Puh^ dy - puh t dy (Al) 


pX 

which is equivalent to Q= q w dx. After rearranging, equation (Al) becomes 

J 0 

4) (A2) 


or, after dividing by 6 / 6, 


q _ r 


T 

x CO 

Poc^oo^cc 5 

rj 
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1 - 




1 , t 2 Jw_u_f 

h t ,ooW 


m: 
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TABLE I.- BOUNDARY CONDITIONS FOR EXPERIMENTAL DATA 

[Pr = 0.725] 


Figure 

Reference 


T t,«’ 

T /T 
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(y - 1)M^ | 

OK 

A W/ A °o 

Poo U oo C p ? oo^t, OO 0 

T* 

w 

N 

00 a 

4uco 

1(a) 

5 

6.78 

639 

4.64 

0.9239 

1.797 

5 

-2.276 

1 



1 

[ 

I 

1 

1 

7 

-2.017 
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t 

T 

j 

♦ 

1 

9 

-1.863 

io>) 



5.01 

399 

4.31 

.4680 

.332 

5 

-.512 

f 



1 

I 

1 

I 

I 
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(a) Moo = 6.78; T w /T„ = 4.64. 
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(b) M„ = 5.01; T w /T„ = 4.31. 

Figure 1.- Experimental temperature-velocity profile data obtained from reference 5 compared with theoretical results calculated by 
using the methods of the present paper, reference 4, and reference 15. 
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(b) Mco = 3.00; T w /Tco = 1.316. 

Figure 2.- Experimental temperature-velocity profile data obtained from reference 6 compared with the present theoretical results. 




u / u / u* 

(C) IVU = 5.01; T w /Too = 4.29. (d) IVU, = 4.93; T w /Tco = 5.42. 

Figure 3.- Experimental temperature-velocity profile data obtained from reference 7 compared with the present theoretical results. 
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(g) M„ = 5.75; T w /T„ = 6.19. 
Figure 3.- Continued. 
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Figure 3.- Concluded. 


(m) Mo, = 8.18; T w /T„ = 6.60. 
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(c) Moo = 5.16; T w /T„ = 3.704. 


(d) Moc = 5.12; T w /T„ = 3.761. 


Figure 4.- Experimental temperature-velocity data obtained from reference 8 compared with theoretical results calculated by using the method 

of the present paper. 
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(k) Mo, = 5.21; T w /T„ = 5.145. 


Figure 4.- Continued. 
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(c) Meo = 1.994; Tw/Too = 1.720. Id) Me = 1.968; T w /T«, - 1.731. 

Figure 5.- Experimental temperature-velocity profile data obtained from reference 9 compared with the present theoretical results 
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(c) Mo. = 3.035; T w /T„, = 2.668. 


(d) NU = 6.500; T w /T„ = 4.100. 


Figure 6.- Experimental temperature-velocity profile data obtained from references 9, 10, and 11 compared with the present theoretical results. 
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(a) Moo = 5.07; T w /T„ = 4.595. 


(b) M„ = 5.20; T w /T=o = 4.615. 



Figure 7.- Experimental temperature-velocity profile data obtained from reference 12 compared with the present theoretical results. 





(c) Moo = 6.02; T w /Tco = 5.16. (d) M„ = 6.02; T w /T«, = 7.38. 

Figure 8.- Experimental temperature-velocity profile data obtained in the Langley 20-inch hypersonic tunnel and experimental data obtained from 

reference 13 compared with the present theoretical results. 
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u/u ro U/U„ 

(C) M« = 0.85; T w /T„ = 1.031. (d) = 0.85; T^Ta, = 1.037. 

Figure 9.- Experimental temperature-velocity profile data obtained from reference 14 compared with the present theoretical results. 
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Figure 9.- Continued. 


(I) Moo = 0.85; T w /Tco = 1.092. 
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(o) M.0 = 1.98; T w /Tco = 1.630. 


Figure 9.- Continued. 
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Figure 9.- Continued. 


(t) Mo, = 2.98; T w /Too = 2.271. 
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Figure 11.- Energy-deficiency distribution across the boundary layer corresponding to the experimental data of figures 4(m) and 4{n). 
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